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Abstract 

In order to check whether odd-numbered Tomonaga-Luttinger ladders are 
dominated by antiferromagnetic correlations associated with gapless spin ex- 
citations, correlation functions of the doped three-chain Hubbard model are 
obtained with the bosonization at the renormalization- group fixed point. The 
correlation of the singlet superconducting pairing across the central and edge 
chains is found to be dominant, reflecting two gapful spin modes, while the 
intra-edge spin density wave correlation, reflecting the gapless mode, is only 
subdominant. This implies that, when there are multiple spin modes, a domi- 
nant superconductivity can arise from the presence of some spin gap(s) despite 

the coexistence of power-law correlated spins. 
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Recently, a wealth of experimental and theoretical results indicate that interacting elec- 
trons on multiple chains, or ladders, are an interesting realm of correlated systems. An 
increasing fascination toward them has been kicked off by an 'even-odd' conjecture by Rice 
et al., who have proposed that the ladder, at half-filling, with even number of chains should 
be a spin liquid reflecting the absence of gapless spin excitations, while odd-numbered chains 
should be antiferromagnetic (AF) reflecting the presence of gapless spin excitations. [I] |3|] 
This is reminiscent of Haldane's conjecture [f||| for the one-dimensional (ID) AF Heisenberg 
model for integer and half-odd-integer spins. 

When the system is doped with carriers, it is usually supposed that an even-numbered 
ladder should exhibit the interchain singlet superconductivity as expected from the persis- 
tent spin gap, while an odd-numbered ladder should have the usual 2kp spin-density wave 
(SDW) reflecting the gapless spin excitations. In ID, an interacting electron system may be 
described by the Tomonaga-Luttinger liquid. || Thus intensive analytical studies have been 
performed by extending the Tomonaga-Luttinger model analysis to the two-chain ladders. 
These analytical calculations support the superconductivity in double chains within the 
perturbational renormalization-group analysis for weak repulsive interactions. [f7[— pTT[] To be 
more precise, the correlation of the interchain pairing is dominant and much stronger than 
that of the subdominant 4k f charge density wave (CDW) in the above calculations. On the 
other hand, numerical calculations performed directly for the two-chain t — J and Hubbard 
models have also been performed, although the phase diagram including the strong-coupling 
regime has not been conclusive. ||12|-[T6[ 



Experimentally, cuprates SrCu 2 3 and Sr 2 Cu 3 5 are investigated as prototypes of two 
and three chain systems, respectively |17|] . The two-chain system indeed shows a spin liquid 



behavior characteristic of a finite spin-correlation length, while the three-chain system shows 
an AF behavior. 

Theoretically, however, whether the 'even-odd' conjecture continues to be valid for triple 
chains remains an open question. In fact, Arrigoni has looked into the triple chains hav- 
ing weak interactions by using the usual perturbational renormalization-group technique to 



conclude that gapless and gapful spin excitations coexist there. |18j Namely, he has actually 



enumerated the numbers of gapless charge and spin modes on the phase diagram spanned by 
the doping level and the interchain electron tunneling strength. He found that, at half-filling, 
one gapless spin mode exists for the interchain hopping comparable with the intrachain hop- 
ping, in agreement with some experimental results and theoretical expectations. Away from 
the half-filling, on the other hand, one gapless spin mode is found to remain at the fixed 
point in the region where the fermi level intersects all the three bands in the noninteracting 
case. ^From this, Arrigoni argues that the spin-spin correlation should decay as a power 
law. 

On the other hand, his result also indicates that two gapful spin modes exist in addition. 
While the existence of a gapful spin mode crudely favors a singlet superconductivity (SS), 
we are in fact faced here with an intriguing problem of what happens when gapless and 
gapful spin modes coexist, since it may well be possible that the presence of gap(s) in some 
out of multiple spin modes may be sufficient for a dominance of superconductivity. This has 
motivated us, in the present work, to actually look at the correlation functions using the 
bosonization method at the fixed point away from half-filling. Since we have the cuprate 
ladder in mind, we concentrate on the open boundary condition (OBC) across the chain, 
where the central chain is inequivalent to the two edge chains. We find that the interchain 
SS pairing between the central and edge chains is the dominant correlation, which is indeed 
realized due to the presence of the two gapful spin modes. On the other hand, the SDW 
correlation, which has a power law for the intra-edge chain reflecting the gapless spin mode, 
is only subdominant. 

We start from the Hamiltonian, 

H = H + H int , (1) 

Ho = ^2 £ka>lrka a irka 
irka 

~ t ^Zi. a ^arku a ^rka + ^rka^rka + h.C.). (2) 
rka 

Here a\ rko . creates an electron with lattice momentum k and spin a on right (r = R) or left 



(r = L) going branch in the 2-th chain (i = a, (3, 7 with (3 being the central one), e k the 
kinetic energy of each chain, and t the interchain hopping. The one-electron part, H , may 
be diagonalized by a linear transformation, 



/ \ 

O'arka 
^(Srka 
darker j 



( 1 V2 
2 2 



JL n L 

V2 U V2 



V 5 



yj 1 
2 2 J 



C 2 rko 



(3) 



resulting in 



Ho = J2i( e k - V2t)c\ rka ci rka + e k c\ rka c 2rku 



rkc 
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(4) 



Hereafter we linearize the band structure around the fermi points as usual and neglect 
the difference in the fermi velocities of three bands, which will be acceptable for the weak- 



hopping case ||19|| . We focus on the case in which all of three bands are away from half-filling. 

The part of the Hamiltonian, H^, that can be diagonalized in the bosonization only 
includes forward-scattering processes in the band picture, and has the form 



Ha = H svin + H ( 



spin T -'-'charge) 
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^charge = f dx [~^{ d xXi+f + K pi (d x X 



(5) 



Here 0j + is the spin phase field of the i-th band, Xi+ is t ne diagonal charge phase field, while 
4>i-(Xi-) is the field dual to 4>i + (xi+), K ai (K pi ) the correlation exponent for the phase 
with v a i(v p i) being their velocities. For the Hubbard type interaction, we have v a i = v F , 
K ai =l for all i's, while v pl = v F , v p2 = v F y/l - 4g 2 , v p3 = v F ^fl - g 2 /4, K pl = 1, K p2 = 
l-2g)/(l + 2g), K p3 = 



l-flf/2)/(l + flf/2), where g = U/2nv F is the Hubbard U 
interaction made dimensionless. 

The diagonalized charge field x%± is linearly related to the initial charge field 9i± of the 
i-th band as 
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where 6i± and 0± are related to the field operator for electrons -?/w as 



T/) i+( _) <T (a;) = ^ )o " exp{±2fcj F a; 



±-{M*) ± + (<&+(» ± &_(*))]}■ 



(7) 



Here rji^'s are majorana fermion operators (or Haldane's U operators) [^0[ which ensure the 
anti-commutation relations between electron operators through the relation, {i]i rcr , i]i' r i cr i} + = 
25u>5 rr i5 acT ', 7)\ ra = rj ira . 

There are still many scattering processes corresponding to the backward scattering and 
pair tunneling scattering processes between two bands, which cannot be treated exactly. 
Arrigoni examined the effect of such scattering processes by diagrammatic perturbational 
renormalization group technique. He found that the backward-scattering interactions within 
the first or the third band turn from positive to negative as the renormalization is performed 
and that the pair tunneling processes between the first and third bands also become relevant. 
At the fixed point the Hamiltonian density, H*, then takes the form, in term of the phase 
variables, 



H* 



^cos(20 1+ (x)) - ^cos(20 3+ (^)) 
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where #&(1), <?&(3) are negative large quantities and 3) is a positive large quantity ]2T 

This indicates that the phase fields <f>i+, 03+, and xi- are long-range ordered and fixed 
at 7r/2, 7r/2, and respectively, which in turn implies that the correlation functions 

that contain <pi_, </> 3 _, and Xi+ fields decay exponentially. The renormalization procedure 
will affect the velocities and the critical exponents for the gapless fields, X2±, X3±> an d <f>2±, 
so that we should end up with renormalized v*'s and K* J s. 



In principle, the numerical values of renormalized t>*'s and K*'s for finite g may be 
obtained from the renormalization equations as has been attempted for a double chain by 



Balents and Fisher p0| , although it would be difficult in practice. However, at least in the 
weak-coupling limit, g — > 0, to which our treatment is meant to fall upon, we will certainly 
have v * ~ vp and K* ~ 1 for gapless modes even after the renormalization procedure. 

Now we are in position to calculate the correlation functions. The two-particle correlation 
functions which include the following two particle operators in the band description show 
power-law decay: 

(1) operators constructed from two operators only in the second band (since the charge and 
spin phases are both gapless, electrons in this band should have the usual Luttinger liquid 
behavior), 

(2) order parameters of singlet superconductivity in the first or third bands, ^i+jq^i-KT)! 

V>3+ta)V>3-i(T) • 

As a result, the order parameters that possess power-law decays should be the following, 
where we also give the exponents: 

(A) The correlations within each of the two edge (a and 7) chains or across the two edge 
chains: 

(a) 2k F CDW, OintraCDW = '^l,( 7 ) +T V ; a(7)-T; C inter CDW = ^afrJ+T ^7( a )-T > 

(b) 2k F SDW, O intra SDW = ^a( 7 )+TV'a(7)-|; O in terSDW = ^a^+^lia)-^ 

(c) singlet pairing (SS), O intraS s = Va(7)+T^a( 7 )-i; O intorS s = WtHT^W-I' 

(d) triplet pairing (TS), O intraT s = Va( 7 )+T^a( 7 )-t; O inteiTS = ^(tHT^tW-T > 

(B) The singlet pairing across the central chain (/3) and an edge chain, O cen traiss — 

^a( 7 )+T V>/9-i- 

In the band picture we can rewright O cen t r aiss as primarily comprising O centra iss ~ 
t/?x+tV ? i— 4. — ^3+1^3-1- We may thus call this paring d- wave-like in a similar sense as in 
the two-chain case, in which a pair is called d-wave when the pairing within the bonding 
band and that within antibonding band enter with opposite signs. PJHj 



Thus the edge-chain SDW correlation has a power-law decay, while the SDW correlation 



within the central chain decays exponentially since it consists of the terms containing </>]_ 
and/or </> 3 _ phases. Although we calculate the case away from half-filling, the SDW corre- 
lation should obviously be more enhanced at half-filling. NMR experiments at half-filling 
fT7| show that the nuclear-spin relaxation rate 1/Ti which is represented by the imaginary 
part of the dynamical susceptibility increase with decreasing temperature for the three-chain 
cuprates in contrast to the two-chain case. This is consistent with the present result, since 
the experiments should detect the total SDW correlation of all the chains. 

Intra- or inter-edge correlation functions have to involve forms bilinear in C2 in eq.(3). 
They are described in terms of the second band 82, which does not contain xi ( ec L-(6)), 
a phase-fixed field. Thus the edge-channel correlations are completely determined by the 
character of the second band (the Luttinger-liquid band), while the other phase fields, being 
gapful, are irrelevant. The final result for the edge-channel correlations at large distances, up 
to 2kp oscillations, is as follows regardless of whether the correlation is intra- or inter-edge: 

(Ocdw(*)O^ dw (0)> ~ x-^ +2K ^~ K ^, 
(Osdw(x)OLw(O)) ~ X ~^ +2K ^-^, 

(Oss(x)Ojs(O)) ~ x^ (7 t +7 t ) " i " ;2 , (9) 

(o TS (x)oyo)> ~ * (3 t +3 & ) ~ 3 & 

By contrast, if we look at the pairing Ocentraiss (#) across the central chain and one of 
the edge chains, this pairing, which circumvents the on-site repulsion and is linked by the 
resonating valence bonding between the neighboring chains, is expected to be stronger than 
other correlations as in the two-chain case. The correlation function for O centra iss ( x ) is indeed 
calculated to be 

(O centralS s(x)OL tralS s(0)) ~ x~^ + ^\ (10) 

In the weak (infinitesimal) interaction limit, all the K*'s tend to unity, where the SS 
exponent becomes as small as 1/2 while the exponents of other correlations tend to 2. Thus, 
at least in this limit, the central SS correlation dominates over the others. The duality 

7 



relation (in which the pairing and density-wave exponents are reciprocal of each other |jTT| ) 
is similar to that in the two-chain case, in which the interchain-SS exponent is 1/2 while the 
exponent of the Ah f CDW is 2. 

In summary, we have studied correlation functions using the bosonization method at 
the renormalization-group fixed point away from half-filling. We found that the dominant 
correlation is the interchain singlet pairing across the central chain and either of the edge 
chains. The key message is that there is an example where the dominance of superconduc- 
tivity only requires the existence of gap(s) in some spin mode, despite the coexistence of a 
power-law spin-spin correlation, when there are multiple modes. It would be interesting to 
further look into how the situation for the single, double, triple, chains crosses over to 
the two-dimensional system. 

We wish to thank E. Arrigoni for sending us his work prior to publication. 
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